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keV sterile neutrinos are important as extensions of the Standard Model of particle physics and 
as serious keV dark matter (DM) candidates (Warm DM) in agreement with both cosmological 
and galaxy observations. We study the possible detection of a keV sterile neutrino through its 
mixing with a light active neutrino in the Rhenium 187 and Tritium beta decays. The electron low 
energy beta spectrum keV < Te [Qp ~ i^s) keV is the region where a sterile neutrino could 
be detected and its mass rus measured, (rris being expected 1 to 10 keV from cosmological and 
galactic observations and theoretical analysis). This region is away from the near end-point Qp 
region suitable to measure the mass of the light active neutrino. In order to analyze the sterile 
neutrino effect, we introduce the dimensionless ratio TZ of the sterile neutrino to the active neutrino 
contribution: TZ exhibits a constant plateau that starts at momentum pe — and sharply drops 
for pe near its maximum value; TZ increases with the mixing angle and decreases for increasing 
sterile neutrino mass rria. The Kurie function is expressed in terms of TZ and reflects its properties. 
For ^*'^Re the p-wave electrons dominate the decay by a factor 10 over the s-wave electrons while 
the Tritium decays exclusively through the s-wave electrons. Two very important experiments are 
running at present: MARE and KATRIN for Rhenium 187 and Tritium beta decays respectively; 
(KATRIN concentrates on the end-point electron region but if appropriately modified could look 
for the keV sterile neutrino as well). The carefully computed decay rates and spectra support that 
the largest realistic ratio TZ is about 10~* for current sterile neutrino models, and the number of 
sterile neutrino events [for typical 10 years of MARE data acquisition] is about 10^ — 10^, which are 
not negligible, [namely, 10^ — 10^ sterile neutrinos in 10^^ — 10^^ events]. A more precise detection 
probability should include a careful analysis of the systematic errors, which is beyond the scope of 
this paper. 
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I. INTRODUCTION 



It is well known that dark matter (DM) is not described by the Standard Model (SM) of particle physics. Many 
extensions can be envisaged to include DM particles, coupled weakly enough to the SM particles to fulfill all particle 
experimental constraints, namely the fact that DM has not been detected so far in any particle physics experiment. 
On the other hand, cosmological and astrophysical constraints such as the ones coming from the dark matter density 
and the galaxy phase space density, or alternatively, the universal galaxy surface density, lead to DM candidates in 
the keV mass scale, namely Warm DM (WDM), refs. A keV mass scale sterile neutrino is the front running 

candidate for WDM. Other WDM candidates in the keV mass scale are: gravitinos, light neutralino and majorons 



Considering the first WDM candidate, sterile neutrinos can be naturally embedded in the SM of particle physics. 
They do not participate in weak interactions, and hence they are singlets of color, weak SU(2) and weak hypercharge. 
One sterile neutrino per lepton family is expected, but only the lightest one (i.e. electron family) has a life time of 
the order of Hubble time and can describe the DM. 

In this work, we consider Rhenium 187 and Tritium beta decay experiments to detect a keV mass sterile neutrino as 
a DM candidate. The left neutrino flavor state Ve (and equivalently for i>e) will be a mixing of two mass eigenstates: 
one light active neutrino mass state (vi) and one keV scale sterile neutrino mass state (I's)- Other neutrino mass 
states will not be taken into account for the time being. The mass mi of the lightest active neutrino state is negligible 
{rni ^ eV ) in comparison with the mass of the keV sterile mass state. The smallness of the mixing angle C makes 
sterile neutrinos difficult to detect. 

Sterile neutrinos in the beta decay of Rhenium 187 are currently searched by the Microcalorimeter Arrays for a 
Rhenium Experiment (MARE) [7]. In this decay the available energy is (5^(^*^Re) ~ 2.47 keV. The beta decay of 
Rhenium 187 into Osmium 187 is a first forbidden unique Gamow- Teller process (5/2+ 1/2^). 

Up to now, the no observation of keV scale sterile neutrinos in the beta decay of Rhenium 187 gave an upper bound 
on the mixing angle ^ < 0.095 for 1 keV steriles Q, which is compatible with the cosmological constraints on the 
mixing angle, ^ ~ 10~^, appropriate to produce enough sterile neutrinos to account for the observed DM. However, 
the amount of the sterile neutrinos that could be produced in the early universe also depends on the production 
mechanism which is model dependent. 

The Karlsruhe Tritium Neutrino Experiment (KATRIN) is currently studying the Tritium beta decay [ll| and, 
if suitably adapted, it could show the presence of a sterile neutrino as well. In this decay the available energy is 
(^^(■^Hi) ~ 18.6 keV. The beta decay of Tritium into Helium 3 is an allowed transition (1/2"*" — > 1/2+) with Fermi 
and Gamow- Teller contributions. Clearly, KATRIN has the potential to detect sterile neutrinos with mass up to 18 



Detection of massive neutrinos by /3-decay has been proposed in ref. [1^ . Other methods proposed to detect sterile 
neutrinos include measurements of the nuclear recoil [is'l and sterile neutrino capture on /3-decaying nuclei [l4| . 




In Sections II and HI we deal with keV dark matter from the cosmological and galactic point of view. WDM (DM 
particle mass between 1 keV and 10 keV, and decoupling temperature Td ^ 100 GeV) produces the observed small 
(galactic) structures, as well as the large scale and cosmological structures, the observed cored density profiles and 
the right surface density value, while GeV WIMPS (m ~ 100 GeV, and ~ 5 GeV, Cold DM) inevitably produce 
a host of non observed small scale structures, non observed cusped profiles and values of the galaxy surface density 
much higher than the observed value. 

In Section IV we analyze the possible detection of a sterile neutrino in the Rhenium beta decay. Contributions for 
the electron s- and p- waves have to be taken into account (9l.[l0j. For the detection of a sterile neutrino, the important 
electron kinetic energy range Te is between and {Qp — Wg), being the mass of the keV sterile neutrino. On the 
contrary, the electron kinetic energy region close to the end-point energy is the one suitable for the detection of 
light active neutrinos. Systematic uncertainties such as Beta Environmental Fine Structure (BEFS) are not considered 
here 0]. In order to analyze the sterile neutrino effect, we introduce the dimensionless ratio Ti. of the sterile neutrino 
to the active neutrino contribution: TZ is largest for Te (or the electronic momentum pe) going to zero. This procedure 
is useful because it allows us to compare two regions of the same spectrum: the one where the keV neutrino imprints 



keV. 
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a glitch on the spectrum and the opposite side (near the end-point) where the active hght neutrino effect shows up. 
For ^^^Re the electrons emitted in the p-wave dominate the decay by a factor 10** over the s-wave. TZ is non-zero in a 
range < Pe < (pe)max, the maximum momentum (pe)max decreases with increasing mass ms- TZ exhibits a constant 
plateau that starts at Pe = and sharply drops for pe near (pe)max- 

TZ increases with the mixing angle and decreases for increasing mass nis . The Kurie function is expressed in terms 
of TZ and the above properties reflect themselves in the obtained Kurie plots. 

In Section V we study the possible detection of a sterile neutrino in the Tritium decay [ll| . Here, only the electron 
s-wave contributes. Analogously to the Rhenium beta decay, the important kinetic energy region for the sterile 
neutrino detection is the low energy range < < (Qp — rus) while for active neutrinos it is the one close to the 
end-point energy Qp. The properties of the ratio TZ and Kurie functions are similar to those obtained in the Rhenium 
case. 

In Section VI we present our conclusions. The low electron energy domain of the beta spectrum is the region where 
a sterile neutrino could be detected and its mass measured, the expected mass being in the keV scale (1 to 10 keV) 
as constrained from cosmological and galactic observations and theoretical analysis. We found that the largest ratio 
TZ for a realistic mixing angle is about lO^'^, [dotted line and small Pe region in Fig. 5]. In order to detect such a 
small sterile neutrino effect, the number of detected events Np must be increased, (for instance by choosing a small 
Qj3 value, as it is the case of the Rhenium, and by increasing the time of data acquisition). We evaluate the sterile 
neutrino probability TZ x Np is about 10^ — 10^, which is not negligible; (for a MARE number of events of 10^'^ — 10^^, 
10 years of data acquisition, 8 arrays and 400 gr of natural Rhenium [3])- That means to find such 10^ — 10^ sterile 
neutrinos in those 10^'^ — lO^'' events. This number of sterile neutrinos is increased to 10^ — 10^ for the MARE option 
of 10^^ events, (10 years of data acquisition, 16 arrays and 3.2 kg of natural Rhenium, [7]). Of course, in order to 
assess a precise detection probability one should include a careful analysis of the systematic errors and instrument 
parameters, but such study is beyond the scope of the present paper. Natural units h = c = 1 are used all over this 
paper. 



II. DARK MATTER 

Although dark matter was noticed seventy-five years ago [l^ [l^ its nature is not yet known. Dark matter (DM) is 
needed to explain the observed structures in the Universe, in particular galaxies. DM particles must be non-relativistic 
by the time of structure formation in order to reproduce the observed small structure at '-^ 2 — 3 kpc. 

The connection between the scale of the formed structure and the mass of the DM particle follows from the value 
of the free-streaming length Ifs [17,]. This is the distance that the DM particles can freely travel. Structures at scales 
smaller than Ifs are erased by free-streaming and hence Ifs provides a lower bound on the size of DM dominated 
structures. WDM particles with mass in the keV scale give Ifs ~ 100 kpc while 100 GeV cold dark matter (CDM) 
particles produce an extremely small Ifs ~ 0.1 pc. A Ifs ^ 100 kpc is in nice agreement with the astronomical 
observations on galaxies [l^ (smaller objects like stars are formed by baryons, not by DM), as well as on cosmological 
scales. 

The GeV CDM Ifs is a. million times smaller and would lead to the existence of a host of CDM smaller scale 
structures till the size of the solar system. No structures of such type have ever been observed. Lighter DM particles 
in the eV scale (hot dark matter, HDM) have a free-streaming length Ifs ~ Mpc and hence would erase all existing 
structures below the Mpc scale in contradiction with all observations. This is why HDM cannot be considered [l9| . 

The reason why CDM does not work is simple: CDM particles in the GeV scale are too slow (too cold) which 
prevent them to erase the small scale structure, while the eV particles (HDM) are excessively fast, which erase all 
structures; in between, WDM keV particles are able to produce the observed structures. 

Astronomical observations strongly indicate that dark matter halos have cored profiles till scales below 1 kpc. On 
the contrary, CDM simulations (particles heavier than 1 GeV) always give cusped profiles. No cusped profiles have 
been ever observed. Linear profiles computed from the Boltzmann-Vlasov equation turn out to be cored for WDM 
and cusped for CDM indicating that WDM does reproduce the astronomical observations Q . 

The surface density in DM dominated galaxies is defined by /io = po rg where pa is the central core density and tq 
is the core radius for DM dominated galaxies, po turns out to be universal, taking the same value up to ±10% for 
galaxies of different sizes, morphologies, Hubble types and luminosities ^20,] . The surface density value predicted by 
CDM simulations is 1000 times larger than the observed value [2l[, while the surface density for keV WDM computed 
from the Boltzmann-Vlasov equation is in full agreement with the observed value of 120 (MeV)^ indicating again that 
WDM does reproduce the astronomical observations ^]. 
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Combining theoretical analysis with the observed values of the dark matter density and the phase space density 
today (density over the cube of the velocity dispersion) of dwarf spheroidal galaxies constrains the mass to be in the 
^ keV range 

Recent radioastronomy observations of velocity widths in galaxies from 21cni HI surveys clearly favours WDM 
over CDM WDM simulations contrasted to astronomical observations suggest a WDM particle mass slightly 

above 1 keV. Constraints from large scale structure give this value too 23]. Recent cosmological WDM N-body 
simulations with keV sterile neutrino WDM clearly show the agreement of the predicted small scale structures with 
the observations, while CDM simulations do not agree with observations at such scales [24i] . 

None of the predictions of CDM simulations at small scales (cusps, substructures, dark disks, ...) have been 
observed. For example, the CDM satellite problem (CDM simulations predict too many satellites in the Milky Way; 
only 1/3 of satellites predicted by CDM simulations around our galaxy are observed) ; the surface density problem 
(the galaxy surface density value for CDM simulations is 1000 larger than observed [l,[2lj); the voids problem and 
the size problem (CDM simulations do not produce big enough galaxies) [2^-[13|- 

Notice that all DM observable effects discussed above only arise from the gravitational behaviour of the DM. Galaxy 
properties are independent of the spin, statistics and non-gravitational couplings of the DM particles, provided that 
their couplings are small enough. 

DM may decouple at or out of thermal equilibrium. The distribution function freezes out at decoupling. Whether 
they decouple at or out of equilibrium depends on the non-gravitational couplings of the DM particle. Normally, 
sterile neutrinos are so weakly coupled that they decouple out of thermal equilibrium. The functional form of the DM 
distribution function depends on the DM particle couplings and is therefore model dependent. 

Sterile neutrinos can decay into an active-like neutrino and a monochromatic X-ray photon with an energy half the 
mass of the sterile neutrino. Observing the X-ray photon provides a way to observe sterile neutrinos in DM halos jl^]. 

WDM keV sterile neutrinos can be copiously produced in the supernovae cores. Supernovae (SN) stringently 
constrain the neutrino mixing angle squared to be < 10~^ for m > 100 keV (in order to avoid excessive energy lost) 
but for smaller masses the SN bound is not so direct. Within the models worked out till now, mixing angles are 
essentially unconstrained by SN in the keV mass range [29j . 

Sterile neutrinos are produced out of thermal equilibrium and their production can be non-resonant (in the absence 
of lepton asymmetries) or resonantly enhanced (if lepton asymmetries are present). keV sterile neutrino WDM in 
minimal extensions of the Standard Model is consistent with Lyman-alpha constraints within a wide range of the 
model parameters. Lyman-alpha observations give a lower bound for the sterile neutrino mass of 4 keV only for sterile 
neutrinos produced in the case of a non- resonant (Dodelson-Widrow) mechanism (30l . [3l| . The Lyman-alpha lower 
bounds for the WDM particle mass are smaller in the Neutrino Minimal Standard Model (sterile neutrinos produced 
by the decay of a heavy neutral scalar) and for fermions in thermal equilibrium. 

Moreover, the number of observed Milky- Way satellites indicates lower bounds between 2 and 13 keV for different 
models of sterile neutrinos. 

In summary, contrary to CDM, WDM essentially works, naturally reproducing the astronomical observations of 
structures over all scales, small as well as large and cosmological scales. The sterile neutrino with mass in the keV scale 
appears as a serious candidate for WDM. Galaxy observations alone cannot determine the DM particle properties other 
than the mass and the decoupling temperature. A direct particle detection is necessary to pinpoint and determine 
which particle candidate describes the dark matter. Beta decay is a promising way to detect DM sterile neutrinos. 



III. DARK MATTER AND KEV STERILE NEUTRINOS 



As is known, DM is not described by the Standard Model of particle physics (SM) . Many extensions of the SM can 
be envisaged to include a DM particle with mass in the keV scale and coupled weakly enough to the Standard Model 
particles to fulfill all particle physics experimental constraints, namely, the fact that the DM has not been detected 
so far in any particle physics experiment. 

Besides sterile neutrinos, possible DM candidates in the keV mass scale are: gravitinos, light neutralino, majoron 

0- ..... . 

As particle physics motivations for sterile neutrinos one can advance that there are both left and right handed 

quarks (with respect to the chirality) while active neutrinos are only left handed. It is thus natural to have right 
handed neutrinos i^fj besides the known left-handed active neutrino. This argument is the so-called 'quark-lepton 
similarity'. 

Sterile neutrinos can be naturally embedded in the SM of particle physics with the symmetry group SU{3)coior ^ 

weak hypercharge- 
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Leptons are color singlets and doublets under weak SU(2) in the SM. Sterile neutrinos vji do not participate to 
weak interactions. Hence, they must be singlets of color, weak SU(2) and weak hypercharge. 

Let us consider a simple embedding of the sterile neutrino in the Standard Model. More elaborated sterile neutrino 
models have been put forward "m] . 

The SM Higgs $ is a SU(2) doublet with a nonzero vacuum expectation value <E>o. This allows a Yukawa- type 
coupling with the left and right handed leptons: 

Lvuk = y VL vr ^0 + h.c. , (1) 

where y = Yukawa coupling, 

<fo = ^ ) , = 174 GeV . 

These terms in the Lagrangian induce a mixing (bilinear) term between vj^ and z/^ which produces transmutations 
i^L ^ i^B.- Mixing and oscillations of particle states are typical of low energy particle physics. Further well known 
examples are: (i) flavor mixing: e-/x neutrino oscillations (explain solar neutrinos), (ii) K'^ — K , — B and — D 
meson oscillations in connection with CP-violation. 

As a consequence of the Lagrangian eq.(IT]), the neutrino mass matrix takes the form: 

where M is the mass term of the right-handed neutrino i^r, rriD = y v and M ^ rriD. 

The masses of the active and sterile neutrinos are given by the seesaw mechanism. The mass eigenvalues in this 
simple model take the form: m\,/M (active neutrino) and M (sterile neutrino), with eigenvectors: 

• active neutrino: VacUve - t^L - ^ I'R ■ 









V ™D 


M I 





• sterile neutrino: v sterile — ^r. + ^rf- ^l, M ^ m|,/M 



M 
n 

2 //\,f oK^n,+ 1 n-5 



Choosing M ^ 1 keV and mo ^ 0.1 eV yields mjj/M about 10^ eV consistent with observations. This corresponds 
to a mixing angle C ^ m,]j/M about 10^'* and would be appropriate to produce enough sterile neutrinos accounting 
for the observed DM. However, notice that the amount of the sterile neutrinos produced in the early universe also 
depends on the production mechanism which is model dependent. The smallness of the mixing angle C makes sterile 
neutrinos difficult to detect. 

Sterile neutrinos in the beta decay of Rhenium 187 are currently searched for by the MARE experiment 0]: 
-'^^^Re — y^^"^ Os-|-e~ -|-z>e . In this decay the available energy: is Q(^®''Re) ~ 2.47 keV. Up to now, the no observation 
of keV sterile neutrinos in the Rhenium 187 beta decay gave the upper bound C < 0.095 for 1 keV sterile neutrinos 
0. This value is largerly compatible with the cosmological constraints on the mixing angle. 

On the other hand, precise measurements of nucleus recoil in Tritium beta decay: ^Hi — >■ '^He 2 + + i> can show 
the presence of a sterile neutrino instead of an ordinary active D in the decay products. In this decay the end-point 
energy is Q('^Hi) — 18.6 keV. The KATRIN experiment is currently studying this beta decay in the end-point region 
P- 

Concerning the number of sterile neutrinos, one sterile neutrino per lepton family is expected. Only the lightest 
one (ie. electron family) has a lifetime ~ Hubble time and can describe the DM. 

In summary, the empty slot of right-handed neutrinos in the Standard Model of particle physics can be filled, in a 
fully consistent way, by keV-scale sterile neutrinos describing the DM. 

IV. THE STERILE NEUTRINO MASS AND THE RHENIUM 187 BETA DECAY 

As a probe to detect possible mixing of keV sterile neutrinos with light active neutrinos, we consider the beta decay 
of Rhenium 187 C^'^Re; Z ==75,A= 187) into Osmium 187 (^^^Os; Z = 76, A = 187), 

i87Re ^ i«^Os + e-+i>e (2) 
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The neutrino flavor eigenstate (and equivalently for v^) can be written as a combination of mass eigenstates of 
light active (subscript i) and heavier mass states. 

|j.e> =^f/edf^^)+5IC^-k.> (3) 

i s 

where the quantities U belong to the unitary leptonic mixing matrix. For the purpose of this paper, we approximate 
this combination as a two mass-eigenstate mixing given by 

|j/e) = cosC \vi) + sinC \vs) (4) 

where is the mixing angle between a light neutrino mass state vi^ and the heavy sterile neutrino mass state Vg. 
Other neutrino mass states will not be taken into account for the time being. 

An effective mass mi can be used for the former combination of light mass active neutrinos, but its value (mi < 
eV) is negligible in comparison with the sterile neutrino mass in the keV scale. 

As for the mixing angle the cosmological constraints based on the observed average DM density suggest Q 

sin^ C - 10"** , C ^ 0.006° . (5) 

However, these constraints on the value of C depend both on the sterile neutrino model and on the sterile neutrino 
production mechanism. Eq. ^ corresponds to currently popular models of DM sterile neutrino [H . 

^^^Re is a long half-life isotope (t\i2 — 4.35 • 10^" years), with ground state spin-parity assignment = 5/2+, 
that has a single ^"-decay branch mode to the ground state l/2~ of ^*^0s, with an endpoint energy Qfj ~ 2.47 keV 
((3/3 ^ Te + m^ + Tu). 

In this transition, the change in total angular momentum is A J = 2 and there is also a change of parity (Att — — ), 
therefore we are dealing at best with a first forbidden Gamow- Teller process. The lepton system (e — v) carries an 
orbital angular momentum L = 1 (first forbidden transition) and a spin S = 1 (unique Gamow- Teller transition) , that 
couple to the total angular momentum J — 2. The two possible angular components of the system, [(?j)e(^j)p] j=2j 
are [(p3/2)e(si/2)p]j=2 and [(si/2)e(P3/2)p] J=2- Therefore, as noted in Q, the total differential decay rate dT/dEe is 
a sum of the two contributions corresponding to the emission of electrons in p-wave and s-wave 

dEe ~ dEe dEe ^ ' 

Following eq. (|3]), we write the theoretical spectral shape of the electron in an (/j)-wave as a sum of the contributions 
from active (light, 1) and sterile (s) neutrinos, 

dTi. dT\. ^ dVl , 

cos\+ sin^C (7) 



where 



dV 



X 



dE,, 



dEe dEe dEe 



= C BRe RRePeP., Ee {Eo - Ee) Fo{Z,Ee) Sl{pe, P.^) e(Eo - Ee - m^) , (8) 



for X = Z,s ;Z stands for the atomic number of the daughter nucleus, FQ{Z,Ee) is the Fermi function and 0{Eo — 
Ee — m^) is the step function; B^e is the dimensionless squared nuclear reduced matrix element (r.m.e), and C 
is a constant, both to be defined later on; Rjie is the nuclear radius 35]. In the above expression, Ee^ Eq and 
Pe = \/ E^ — TOg are the total energy, maximum total energy and momentum of the emitted electron respectively, and 
Pi, = y/ {Eq — Ee)"^ + m^ is the momentum of the emitted neutrino. 

Being Qp the end-point energy, we have Eq = me + Qp for zero neutrino mass, and the kinematical ranges of Ee^ Pe 
and are: 



me < Ee <me + Qp; <Pe< \ Ql + '2 meQp; < Pu < Qp- (9) 



The shape factor Si{pe,Pi^) appears in forbidden decays and for the case of interest here, a first forbidden decay, 
takes the value I = for the s-wave and / = 1 for the p-wave electrons. We have 
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1 1 (ZE^ 

So{Pe,Pu) = izpI and Si{pe,Py) ^ izpI ^ J". ■ (10) 

The relativistic Fermi functions Fo{Z, Ee) and Fi{Z, Ee) account for the Coulomb interaction between the residual 
nucleus (Z — 76 in our case) and the emitted electron in the s and p-waves respectively. They are defined as 



Fk-i = 



r(2fc + i) "'^ 



_rwr(i + 2..)J (^P. «)''^-" I rb. + (n) 

and depend on the strength of the Coulomb interaction, given by the fine structure constant a ~ 1/137.03, through 

Ee 



7fc Jk"^ - iaZy and z^aZ — , (12) 

Pe 

k = 1, 2 in our case. We note that the Fermi functions in eg. pTj) satisfy Fi;_i{Z — > 0,i?e) — > 1 for a Z — > and for 
any fc > 1. 

The constant factor C, in eq. ([8]), is given by 

C^^^Flk±^2xlQ-^^ {kc\y\ (13) 

2 TT'' 

where Gp is the Fermi constant, Vud the element of the Cabibbo-Kobayashi-Maskawa matrix (| Vud \— 0.97), and 
~ 1 is the vector strength of the charged weak interaction. 

The dimensionless squared nuclear reduced matrix element (r. m. e.) Bue is defined as 

Si^e - ^ |( i«^0s(l/2-) II ^ r+K- ® r^^ \\ i«^Re(5/2+) ) p (14) 

Re j^i 

(6 = 2 J + 1 as J = 5/2) with [(jj (g) rj]2 a rank-2 tensor given by: 

[CT(g)r]2^= ^ {lplq\2p) (Tprq (15) 

p,g=±l,0 

where /i = ±2, ±1,0; tfj, denote nucleon spin and space coordinates, respectively, and t'^ is the isospin operator 
that transforms one neutron into one proton. The reduced matrix element can be computed directly from the 
experimental ^^^Re mean-life r — ti/2/\n2, [eq. (|8|)] 

B^l ^tC / Pe p. {Eo - E,) Fo{Z, E,) S{pe,p,) dE, , (16) 

audit takes the value i3_Re — 3.6 x 10^'' for a value of the nuclear radius Rue approximated as RRe — 1.2 x (187)^/'^fm ~ 
6.86 fm. 

In Fig. [T] we represent on the left dT\_^^^l dE^ and dr^^^^/dEe, and on the right dVp^^^/dEe and dTp^^^/dEf,; all 
these decay rates are normalized to one. We plot the sterile neutrino contribution, for s and p-waves outgoing single 
electron, for selected values of the sterile neutrino mass, rus =1, 1.5 and 2 keV, compared to the light neutrino case 
with mi = 0. Solid line corresponds to mi = 0; and dashed, dash-dotted, dotted lines correspond to the sterile masses 
respectively. 

In Fig. [21 we represent the non normalized decay rates: on the left dr[_^^^/dEe and dTl^^^/dEe, and on the right 

dTp^^^/dEe and dVp^^^/ dE^. Analogously to Fig. [T] the same neutrino masses are considered. As seen in the plots, 

the maximum differential decay rate for s-wave electrons is of the order of 10^, whereas for p-wave electrons it is 10^"'^. 
This dominance by four orders of magnitude of the p-wave is true for both light and sterile neutrino emission, and 
it has been noticed both theoretically [9] and experimentally [l0|. This is why the spectral shape of beta decay is 
dictated by the curves shown in the r.h.s panel of Fig. [51 

The effect of the sterile neutrino emission on the electron spectral shape is qualitatively represented in Fig. [31 
through the differential decay rate dT/dEe for mg = 1 keV and the dominant contribution cos^ C {dV^ / dE^). To make 
the effect more dramatic a large unrealistic mixing angle sin^ C, = 0.4 is chosen. The two curves start to deviate at the 
step point — Qp — = 1.47 keV, where the sterile neutrino starts to contribute, and the difference grows when 
Te goes to zero. 
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FIG. 1: Contributions of s-wave (left) and p-wave (right) electrons to the normalized differential decay rate of the process ^*^^Re 
to ^*^Os against the electron kinetic energy Ee — irie. Selected values of the sterile neutrino mass are used, rus = 1, 1.5, 2 keV 
(dashed, dashed-dotted and dotted lines, respectively), compared to the light neutrino case m; = (solid line). 

A. The ratio of the sterile neutrino to the active neutrino contribution and the Kurie plot 

In order to analyze the possible effect of a sterile neutrino, we introduce the dimensionless function 

that weights the relative contribution of the sterile neutrino to the light neutrino, times the tangent square of the 
mixing angle. The function TZ is largest for pe (or Te) going to zero. This procedure is useful because we are 
comparing two regions of the same spectrum: the region where {Ef. — me) < (Qp — mg) and the emitted neutrino has 
enough energy such that the sterile neutrino (m^ ^ keV) imprints an effect on the spectrum, and the region where 
(-Ee — nrie) > (Q/3 — TOs) ^ud the sterile neutrino effect does not show up. Clearly, there is a step in the spectrum for 
Ee — me = Qi3 — TTis which could be observed if the experimental relative error in this energy region is lower than the 
height of the step. 

The ratio TZ [eq. (I17p ] is shown in Fig. 2] as a function of the electron momentum pe for a fixed (unrealistic) mixing 
angle, ^ = 0.1°, and for different values of the neutrino masses: mg = 0, 1, 1.5 and 2 keV, corresponding to the solid, 
dashed, dash-dotted and dotted lines respectively. As can be seen in this figure, the ratio is different from zero in the 
rangeO <Pe < (Pe)max, with (pe)max = [(Q/j - TTis) (Q^ - TTis + 2me)] . For example: for = 2 keV, {pe)max — 21.7 
keV; and for ms = 1 keV, {pe)max — 38.4 keV. Notice that when increases TZ decreases. 

Similarly, in Fig. [5] we show the ratio 72. as a function of the electron momentum pe but for a fixed sterile neutrino 
mass of ms — 1 keV and different light-sterile mixing angles C : C, = 0.1°, 0.05°, 0.01°, 0.005°. Figure [S] shows not 
only the increase of TZ with increasing mixing angle for a fixed value of ms, but also shows the fact that for fixed 
values of and C, the ratio TZ is almost constant in the region < Pe < (Pe)max- 

From eq. ([6]) and eq. ((T)), we write 
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FIG. 2: Contributions of s-wave and p-wave electrons as in Fig. [l]of the process ^^'^Re to ^*^Os but for non-normalized decay 
rates in units of the constant Cije = C Bue Rrs in keV''. The maximum differential decay rate for s-wave electrons is of the 
order of 10^, whereas for p-wave electrons it is 10^^. The p-wave dominates by four orders of magnitude for both light and 
sterile neutrino emission and so the spectral shape of beta decay is dictated by the r.h.s. panel. 



being dV/dEe the total differential decay rate where neutrino mixing is present {( 0), and dT^ /dEe the differential 
decay rate for the light neutrino (no mixing: = 0). For small mixing angle C, the differential decay rate dT/dE^ 
[eq. ([18])] normalized to dV'- /dE^ is 



dr/dEe , , 

which shows that for small mixing angle the ratio between the differential decay rate [with mixing and without 
mixing] is given by 1 + 7?.. This ratio will be largest for pe or Tg going to zero. 

We want to emphasize that the energy region suitable for creation and detection of the keV sterile (non relativistic) 
neutrino corresponds to low Pe or Tg. On the contrary, active neutrinos information should be obtained from the 
region of Tg close to the end-point energy Q^. 

From eqs. (O, (|H]) and (fTUl) . we can write for the function the explicit expression, 



1 



Vi F,{Z,E,) 



-77 am ^ \ < Fo{Z,E^) ^ 

Pi^l ^ , Pe_ J'U^, -C/ej 

pI^ Fo{Z,Ee) 

In order to analyze the ratio Fi{Z, E^)/ Fo{Z, Ee), it is worth to define Fk-i{Z,Ef.) as 



Fk-i{Z,E,) = Ck-idk^i [ — ] [ — aZ] (21) 

Pe J \me 



10 




1 1.5 
E -m [keV] 



2.5 



-1 O '7 1 ft T 

FIG. 3; Qualitative graph showing the differential decay rate of the process Re to Os vs. the electron kinetic energy 
considering light-sterile mixing for rria = 1 keV and unrealistically large ^ (solid line) in order to visualize the sterile neutrino 
effect. The dashed line corresponding to cos^ {dV' /dEe) without sterile component is shown for comparison. 



where 



Ck-i = 2 7r(2mei?) 



2(7fc-fc) 



r(2fc + 1) 



r(fc) r(i + 2 7fc) 



dk-1 = t;- [ — ^ 

2tT \pe 



1-2-tk 



TT a. Z 



(22) 



7fc is defined by eq. (|12l) and dk-i{a Z Ee/pe) ^ 1 for (for all k). The above definitions yield again 

Fk-i 1 for a Z — > (for all k). In this respect, eqs. (|2T|) and (|22]) differ from references @ and [3^|3J|. Finally, 
the ratio in the shape factor for Z = 1 [eq. pU]) ] is given by 



plFi{Z,Ee) _ Ci di 



FoiZ,Ee) 



= — — [— a Z 

Co do \m, 



E, 



2 (72-71) 



(23) 



and for pg — >■ 



plF,iZ,E,) p^^o 
FoiZ,Ee) 

[di/do 1 for pe 0]. From eq. ([20]) and eq. ([24]), we have 



0.4 , 



(24) 



TZ ^ tan 4 



(25) 



Forpe 0, the maximum neutrino momenta are: {pus)ma.x ^ IA5 keV, 2.26 kcV for nis ~ 2 keV, 1 keV respectively; 
and {pu,)max = 2.47 keV for m; = keV [being {pe)min — 0]. 
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FIG. 4: Ratio TZ [eq. (I17|l ] of the sterile neutrino to the active neutrino contributions of the process ^*^Re to ^*'^0s vs. the 
electron momentum for a fixed (unrealistic) mixing angle C, = 0.1° and different sterile neutrino masses, iris ~ 0, 1, 1.5 and 2 keV 
(dashed, dashed-dotted and dotted lines, respectively). TZ increases with decreasing ms. TZ is nonzero in a range < pe < Pmax 
and Pmax decreases as ms increases. 



We have shown the relevance of the function TZ in order to analyze the effect of a sterile neutrino. One can also 
study the effect of a sterile neutrino through the difference between the decay rate with mixing (^ ^ 0) and the 
reference case without mixing {( — 0). This difference is very small, since the mixing is in any case small, as it is 
expressed in the following ratio: 



TZ* = 



' dr ' 




' dV ' 




dEe. 




















C=o 





= -1 



dT'/dEe \ ,.^2 ^ ^ dT/dEe 
dV^/dEj ^^"^ dTi/dEe 



- 1 



which can be written as well as a function of TZ 

TZ* = ~sin2c + 7^cos2C 

In Fig. [6] we plot the quantity TZ* for a fixed sterile mass mg — 1 keV and for different mixing angles. 
The Kurie function is defined as 



(26) 



(27) 



Kiy) = 



dV/dEe 



, y = Eo-E, = Q-T,>Q 



Considering the mixing between the light and sterile neutrinos, K{y) can be written as 



(28) 



K{y) = JKfiy) cos^ C + K^iy) sin' ( 



(29) 
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FIG. 5: The ratio TZ as in Fig. 2] but for a fixed sterile neutrino mass rris = 1 keV and different mixing angles, 
C = 0.1°, 0.05°, 0.01°, 0.005° (solid, dashed, dashed-dotted and dotted lines, respectively) of the process ^^'^Re to ^*^Os 
. TZ is almost constant in the range < Pe < Pmax and increases with ^. 



where 

For C = (no sterile-light neutrino mixing), and due to the introduction of FQ(Z,Ee) hi the denominator, K(y) vs. y 
is a straight line for mi ~ 0. K{y) can be written as well in terms of TZ eq. ([T7|. 

K{y) = Ki VTTn COS C . (31) 
Finally, in Fig. [7] we present the Kurie plot for Kg considering several sterile neutrino masses m^, and Ki (solid line). 



V. THE STERILE NEUTRINO MASS AND THE TRITIUM BETA DECAY 



Let us now consider the beta decay of Tritium -^H (Z = 1; A = 3) 

^ ^He + e- +De (32) 

as a probe to detect a possible mixing of keV sterile neutrinos with active neutrinos. "^H is a hydrogen isotope going 
to the helium isotope '^He, with a half-life ti/2 — 12.33 years, end-point energy ~ 18.59 keV, and a spin-parity 
transition 1/2+ ^1/2+. 

For the Tritium decay case "^H, there is no change in angular momentum and parity corresponding to an allowed 
transition {L = 0) with Fermi {S = 0) and Gamow- Teller {S = 1) components. Therefore, the single electron is 
emitted in s-wave and the differential decay rate is simply 

d^_d£s^ (33) 
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FIG. 6: Same as in Fig. [5] but for TV [eq. (|26|) or eq. (|27|) ] in units of 10 ^ for fixed ms and different mixing angles of the 
process ^^■^Re to ^*^0s. The small Pe region is always the best for the sterile neutrino detection. 



Similarly to eq. ^ and eq. ([8]) we have 
and 

— §^ = CBt Pe P., (i^o - Se) Fo{Z , Ee) e{Eo - Ee ^ m^) , x = l,s , (35) 

as the shape factor S{pe,Pi^) — 1 for allowed decays. Here, the squared reduced matrix element (r.m.e) for the allowed 
decay of Tritium (T), will be 

Bt = Bf + Bgt (36) 
where Bp and Bqt are the Fermi and Gamow- Teller decay strengths given by: 

A— 3 

Bf^\ |( =^He(l/2+) II ^ r+ || ^^{1/2+) )\^ (37) 

and 

Bgt - ^ |('He(l/2+) || X^^/'^l II 'H(l/2+))p (38) 

where 17^1 = ca/ cv — 1.26 is the axial-to- vector strength ratio of the charged weak interaction. From the experimental 
half-life of Tritium one obtains a value of the squared r.m.e. Bt — 5.61. 
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FIG. 7; Kurie plot Ks of the process Re to Os for sterile neutrino masses of ms = 1, 1.5 and 2 keV as well as Ki for a 
light neutrino mass mi — 0. 



In Fig. [8] we plot the differential decay rates of the process ■^H to ^He, normalized to 1. We consider the emission of 
a light neutrino, mi — (solid line); and a sterile neutrino with masses = 1, 1.5 and 2 keV (dashed, dashed-dotted 
and dotted lines respectively). 

In Fig. iniwe plot the ratio TZ vs. pe for a fixed sterile neutrino mass ( rn^ — 1 keV) and different mixing angles. 



VI. CONCLUSIONS 



The detection of sterile neutrinos is not only important from the point of view of particle physics for the extension 
of the SM, but also from the point of view of cosmology and astrophysics as a serious candidate for dark matter in 
the keV mass range. With the relevance of the search and detection of keV scale dark matter candidates in mind, 
we have studied Rhenium 187 and Tritium beta decays. The low electron energy domain of the beta spectrum is the 
region where a sterile neutrino could be detected and its mass measured, the expected mass being in the keV scale 
(1 to 10 keV) as constrained from cosmological and galactic observations and theoretical analysis. The electron low 
energy region suitable to detect the sterile neutrino, keV J; 71; < {Qp ~ ms) keV , is away from the near end-point 
energy range suitable for the detection of the active neutrino mass. 

Two very important experiments are running at present: MARE and KATRIN dealing with Rhenium 187 and 
Tritium beta decays respectively. The MARE experiment will provide the entire shape of the electron differential 
decay rate, giving data for both the sterile and the active neutrino detection regions. KATRIN so far concentrates 
on the region near the end-point of the electron spectrum but it would be extremely interesting to also look for data 
for the sterile neutrino in the mass keV region. In this paper, we have carried out the study for the detection of a 
sterile neutrino in the expected keV mass range considering different mixing angles, according to astronomical and 
cosmological observations and experiments. 

In order to detect the small relative deviation in the experimental spectrum due to the sterile neutrino mixing, the 
experimental relative error (inversely proportional to the square root of the number of measured events, e ^ No ) 
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FIG. 8: Differential decay rate for the Tritium process '^H going to '^He with the emission of active (m; = 0) or sterile {iris = 
1, 1.5, 2 keV) neutrinos, in units of the constant Ct = C Bt in keV*; the vertical axis is in logarithmic scale. 



must be as small as possible. In order to reduce the experimental relative error, the number of detected events must 
be increased, for instance by choosing a beta decay with a small Qp value, as it is the case of the Rhenium, and by 
increasing the time of data acquisition. For MARE, the typical number of events is 10^'^ — 10^*, (for 10 years of data 
acquisition, 8 arrays and 400 gr of natural Rhenium [7]). We found that the largest ratio TZ for a realistic mixing 
angle is about 10~^, [dotted line and small region in our Fig. 5 ], therefore the sterile neutrino probability TZ x Np 
is about 10^ — 10^, which is not negligible. That means to find such 10^ — 10^ sterile neutrinos in those 10^'^ — 10^''' 
events. This figure is increased in an order of magnitude for the MARE option of 10^^ events, (for 10 years of data 
acquisition, 16 arrays and 3.2 kg. of natural Rhenium, Q). A simple estimate requires the Poisson error e to be 
smaller than the ratio TZ. Namely, Np > l/TZ? ^ 10^"* — 10^^. Of course, in order to assess a precise prediction of the 
detection probability one should include a careful analysis of the systematics errors and instrument parameters, but 
such study is beyond the scope of the present paper. 
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